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Abstract 

Using an isospin- and momentum-dependent transport model, we study the effects due to the 
momentum dependence of isoscalar nuclear potential as well as that of symmetry potential on 
two-nucleon correlation functions and light cluster production in intermediate energy heavy-ion 
collisions induced by neutron-rich nuclei. It is found that both observables are affected significantly 
by the momentum dependence of nuclear potential, leading to a reduction of their sensitivity to 
the stiffness of nuclear symmetry energy. However, the t/ 3 He ratio remains a sensitive probe of 
the density dependence of nuclear symmetry energy. 
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I. INTRODUCTION 



The equation of state (EOS) of an asymmetric nuclear matter with unequal numbers of 
protons and neutrons depends crucially on the nuclear symmetry energy Although the nu- 
clear symmetry energy at normal nuclear matter density is known to be around 30 MeV from 
the empirical liquid-drop mass formula , its values at other densities are poorly known. 
Studies based on various theoretical models also give widely different predictions [3J]. Lack of 
;his knowledge has hampered our understanding of both the structure of radioactive nuclei 
I fl fl 3 and mmy i m po rt an t 1S8 ne 8 in nnc.ear a S ,oph y8ic8 RflQ, sncn as the nnc,e- 
osynthesis during the pre-supernova evolution of massive stars and the properties of neutron 
stars. However, recent advance in radioactive nuclear beam facilities provides a unique op- 
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portunity to study the density dependence of nuclear symmetry energy 
Theoretical studies have already shown that in heavy-ion collisions induced by neutron-rich 
nuclei, the effect of nuclear symmetry energy can be studied via the ratio of pre-equilibrium 
neutrons and protons [l^, the isospin fractionation [lfiHl?! lsl llil I20I] . 
tifragmentation l2ll 

xansverse flow [2 3! 24^ , the ir~ to tt + ratio [jj , and isospin diffusion in heavy-ion collisions 
2^1 I^tJ . Also, it was found that the correlation function between nucleon pairs with high 



the isoscaling in mul- 
the proton differential elliptic flow |22j , the neutron-proton differential 



jy neutron-rich 
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total momentum and the isobaric ratio t/ 3 He in heavy- ion collisions induced 
nuclei are sensitive to the density dependence of nuclear symmetry energy 
implying that the space-time properties of neutron and proton emission sources are affected 
by the nuclear symmetry energy. In these studies, the momentum dependence of nuclear 
mean-field potential, especially its isovector part (the symmetry potential), was not taken 
into account. The momentum dependence of nuclear isoscalar potential is well-known, and 
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a& lag, 

the experimen- 



ts effect in heavy-ion collisions is large 32]. As shown in Refs. 
tal data on the nucleon directed flow in intermediate energy heavy-ion collisions, which were 
previously explained by a stiff nuclear equation of state with compressibility of 380 MeV 
when the momentum dependence of nuclear isoscalar potential was not taken into account, 
are actually consistent with a soft nuclear equation of state with compressibility of 200 MeV 
after including the momentum dependence. The momentum dependence of nuclear isoscalar 
potential was also found to affect the space-time properties of nucleon emission source [37 1. 
The momentum dependence of nuclear symmetry potential is, on the other hand, poorly 
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known, 
studied 



ts effect on some of the isospin observables mentioned above has only recently been 
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and was found to be significant. In this paper, we shall study how the 
momentum dependence of both isoscalar nuclear potential and nuclear symmetry potential 
affect two-nucleon correlation functions and light cluster production in intermediate energy 
heavy-ion collisions induced by neutron-rich nuclei. Our results show that the momentum 
dependence of nuclear potential reduces the sensitivity of two-nucleon correlation functions 
and the light cluster yield on the stiffness of nuclear symmetry energy. However, the t/ 3 He 
ratio is still found to depend on the stiffness of nuclear symmetry energy even with the 
inclusion of its momentum dependence. 

This paper is organized as follows. In Section |HJ we discuss the momentum dependence 
of nuclear isoscalar potential and nuclear symmetry potential. In Section 11111 we present 
results from the isospin- and momentum-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) 
transport model on two-nucleon correlation functions in heavy-ion collisions induced by 
neutron-rich nuclei at intermediate energies. Results on light cluster production and the 
t/ 3 He ratio based on the nucleon coalescence model are given in Section HVl Finally, we 
conclude with a summary in Section El 

II. MOMENTUM DEPENDENCE OF NUCLEAR MEAN-FIELD POTENTIAL 

The energy per nucleon in an asymmetric nuclear matter is usually expressed as 

E(p, 6) = E(p, 6 = 0) + E sym ( P )5 2 + 0(5 4 ), (1) 

where p = p n + p p is the baryon density with p n and p p denoting the neutron and proton 
densities, respectively; 6 = (p n — p P )/(p p + p n ) is the isospin asymmetry; and E(p, 5 = 0) is 
the energy per particle in a symmetric nuclear matter, while E sjm (p) is the nuclear symmetry 
energy. Studies based on various many-body theories using non-local interactions have shown 
that the momentum dependence of nuclear single-particle potential is different for neutrons 
and protons in asymmetric nuclear matter, see e.g., Ref. |9j for a review. Using the Gogny 
effective interaction (MDI), the nucleon single-particle potential was recently determined in 
the mean-field approximation by fitting the saturation properties of nuclear matter at zero 
temperature with compressibility Kq = 211 MeV and symmetry energy 30 MeV j^]. The 
resulting potential U(p, 5, p, r) for a nucleon with isospin r (1/2 for protons and -1/2 for 
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neutrons) and momentum p in an asymmetric nuclear matter with isospin asymmetry S and 
density p, can be parameterized as 

U(p,6,p,T) = A u ^ + A^ + b(P-\ (l-^ 2 )-8rx B /"'V ' 
Po Po VPo/ o-+l Po 

+ Po 7 p i + (p-p') 2 /a 2+ po y dp i + (p-p') 2 /A 2 ' l j 

where r^r' and / r (r, p) denotes the phase-space distribution function at coordinate r and 
momentum p. The parameter x in Eq.(j2j) is introduced to reflect the uncertainty of our 
knowledge on the density dependence of nuclear symmetry energy, and it has a value of 1 
in the Cogny interaction [4]]. For the momentum- independent part of the potential, given 
by the first line of Eq.(J2J), the value for a is 3/4 and for B is 106.35 MeV, while those for Ai 
and A u depend on the value of x. For x = 0, we have Ai = —120.57 MeV and A u = —95.98 
MeV. For other values of x, they are 

2,B 2,B 

Ai{x) = MO) + — — x, A u (x) = A u (0) —x, (3) 

a + 1 o~ + 1 

with Ai{0) and A u {0) denoting the values for x = 0. 

The terms with parameters C un nk e = —103.4 MeV and Cuke = —11.7 MeV with A = 
l.Opp, where p° F denotes the Fermi momentum at normal nuclear matter density, in the 
second line of Eq. (J2J) describes the momentum dependence of not only the nuclear isoscalar 
potential but also the nuclear symmetry potential as a nucleon with isospin r interacts 
differently with unlike and like nucleons in the background fields. 

With different x values, one can obtain different density dependence for the nuclear 
symmetry energy without changing other properties of the asymmetric nuclear matter. As 
seen in Fig. where the density dependence of nuclear symmetry energy from the MDI 
interaction with different x values are shown, the symmetry energy becomes stiffer with 
decreasing x but its value at normal nuclear density remains the same. To explore the large 
range of E sym (p) predicted by various many-body theories (5, 42, 43, 44 1, we consider in the 
present study two values of x = 1 and x = —2. Since the nuclear symmetry energy given by 
x = 1 has a weaker dependence on density than the one given by x = —2, we call it a soft 
symmetry energy while the one corresponding to x = —2 is called a hard symmetry energy. 
We note that the soft symmetry energy with x — 1 is the same as that given by the default 
Cogny interaction. 
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FIG. 1: (Color online) Density dependence of nuclear symmetry energy from the MDI interaction 
with different x values. 

The difference in the density dependence of soft and stiff symmetry energies is due to the 
different density dependence in the momentum-independent part of the potential from the 
MDI interaction as shown in Fig. |21 where we give the density dependence of total symmetry 
energy for these two cases together with that due to contributions from the kinetic energy 
part as well as the momentum-independent and momentum-dependent potential energy 
parts. The total potential energy contribution to the symmetry energy shown in FigJ2]can 
be parameterized as 

E^Kp) = 3.08 + 39.6m - 29. 2u 2 + 5.68u 3 - 0.523u 4 (MeV) , (4) 

for the soft symmetry energy, i.e., the MDI with x = 1 and 

E p s f m (p) = -1.83 - 5.45m + 30.34m 2 - 5.04m 3 + 0.45m 4 (MeV) (5) 

for the hard symmetry energy, i.e., the MDI with x = —2. In the above, u = p/po is the 
reduced nucleon density. 

To study the effects due to the momentum dependence of nuclear symmetry energy, we 
also consider a nucleon potential U noms (p, 6, p, r) = U (p, p) + U sym (p, 5, r) with the isoscalar 
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FIG. 2: (Color online) Density dependence of total symmetry energy together with contributions 
from the kinetic energy part as well as the momentum-independent and momentum-dependent 
potential energy parts from the MDI interaction with (a) x = 1 and (b) x = —2. 



part taken from the original MDYI interaction 15411 . i.e., 

f/o(p,p) = -H0.44 M + 140.9n 1 - 24 -— / rfV S ^ (6) 

po J 1 + (p - p')7(1.58^) 2 

which has a compressibility Kq = 215 MeV and is almost the same as the momentum- 
dependent isoscalar potential given by the MDI interaction, i.e., Eq. (J2J). For the momentum- 
independent symmetry potential U sym (p, 5, r), it is obtained from U sym (p, 6, r) = dW sym /dp T 
using the isosp in-dependent part of the potential energy density W sym = E^ y ^ n (p) ■ p-5 2 where 
Esymip) IS given by Eqs (J1J and (0) from the MDI interaction. 

To study the effects due to the momentum dependence of isoscalar nuclear potential, 
we also consider the usual momentum-independent soft nuclear isoscalar potential with 
K = 200 MeV (SBKD), firstly introduced by Bertsch, Kruse and Das Guptaj3|, i.e., 

U(p) = -356 u + 303 u 7/6 . (7) 

The density dependence of neutron and proton single-particle potentials at isospin asym- 
metry 5 = 0.2 and at momenta k = 1 fm _1 and k = 3 fm _1 are shown in Fig. El for the 
SBKD, MDYI, and MDI interactions. It is seen that neutrons generally have a stronger 




FIG. 3: (Color online) Density dependence of neutron and proton single-particle potentials at 
different momenta in asymmetric nuclear matter with isospin asymmetry 5 = 0.2 for different 
nucleon effective interactions. 

repulsive potential than protons. Their difference further depends on the nuclear symmetry 
energy, with the soft one giving a larger difference at low densities while the hard one giving 
a larger difference at high densities. Also, with momentum-dependent MDYI and MDI in- 
teractions, the nuclear potential is more repulsive for high momentum nucleons. Moreover, 
the momentum dependence of nuclear symmetry potential from the MDI interaction gen- 
erally reduces the difference between the neutron and proton potentials, especially at high 
momenta. 
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FIG. 4: (Color online) Momentum dependence of neutron and proton single-particle potentials in 
asymmetric nuclear matter with isospin asymmetry <5 = 0.2 and at different densities for different 
nucleon effective interactions. 
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To see more clearly the momentum dependence of nucleon mean-field potential, we show 
in Fig. |3] the momentum dependence of neutron and proton single-particle potentials in 
asymmetric nuclear matter with isospin asymmetry 5 = 0.2 and at different densities for dif- 
ferent nucleon effective interactions. Although the MDI interaction is seen to give a similar 
momentum dependence in the nucleon potential as the MDYI interaction, the difference be- 
tween the neutron and proton potentials, is reduced, especially at higher momenta, due to its 
momentum dependence in the nuclear symmetry potential. Therefore, the MDI interaction 
gives a weaker symmetry potential than that from the MDYI interaction. 
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FIG. 5: (Color online) (a) Density and (b) isospin asymmetry dependence of nucleon effective 
mass in asymmetric nuclear matter with isospin asymmetry 5 = 0.2 and at normal nuclear matter 
density for the MDYI and MDI interactions. 

The nucleon effective mass is an important quantity to characterize the momentum de- 
pendence of nuclear mean-field potential. From the single-particle potential, we can obtain 
the nucleon effective mass from 



m; 



1 + 



m T dU 



(8) 



P d PS P =Pf(T) 

In Fig. El we show both the density (a) and isospin asymmetry (b) dependence of nucleon 
effective mass in asymmetric nuclear matter with isospin asymmetry 5 = 0.2 and at normal 
nuclear matter density for the MDYI and MDI interactions. It is seen that a neutron has 
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a larger effective mass than a proton in neutron-rich nuclear matter and both depends 
strongly on the isospin asymmetry of nuclear matter. The small difference between the 
neutron and proton effective masses for the MDYI interaction is due to the different neutron 
and proton Fermi energies. In addition, the MDYI interaction shows a slightly stronger 
density dependence for the nucleon effective mass than the MDI interaction. 

Above discussions indicate that the nucleon symmetry potential is reduced after including 
its momentum dependence and the effect gets larger with increasing nucleon momentum. 
The difference between symmetry potentials with and without momentum dependence is 
thus larger for nucleons with higher momenta. In the following, the three potentials ob- 
tained with the SBKD, MDYI, and MDI interactions will be used together with the isospin- 
dependent experimental NN cross sections in the IBUU model 



15 



to study the 



effects due to the momentum dependence of isoscalar and isovector (symmetry) potentials 
on two-nucleon correlation functions and light cluster production in intermediate energy 
heavy- ion collisions induced by neutron-rich nuclei. In previous studies 3, \lL Q ^ has 
been shown that these observables are sensitive to the symmetry potential but not to the 
isoscalar potential and the in-medium nucleon-nucleon cross sections. 



III. TWO-NUCLEON CORRELATION FUNCTIONS 



The space-time properties of nucleon emission source, which are important for under- 
standing the reaction dynamics of heavy-ion collisions, can be extracted from two-particle 



correlation functions; see, e.g., Refs. |4( 
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for reviews. In most studies, only 
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561 ] . Recently, data 



the two-proton correlation function is studied 
on two-neutron and neutron-proton correlation functions have also become available. The 
neutron-proton correlation function is especially useful as it is free of correlations due to 
wave-function anti-symmetrization and Coulomb interactions. Indeed, Ghetti et al. have 
deduced from measured neutron-proton correlation function the emission sequence of neu- 



58, 



and have also 



trons and protons in intermediate energy heavy-ion collisions J5j 
studied the isospin effects on two-nucleon correlation functions 601. 

In standard Koonin-Pratt formalism 0, 0, Q, the two-particle correlation function 
is obtained by convoluting the emission function g(p,x), i.e., the probability for emitting 
a particle with momentum p from the space-time point x = (r, t), with the relative wave 
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function of the two particles, i.e., 

r(p _ Jrf 4 x 1 ^(P/2,x 1 )^(P/2,x 2 ) |0(g,r)| 2 
1 ' qJ fd*x 1 g(P/2,x 1 )Jd*x 2 g(I>/2,x 2 ) ' U 

In the above, P(= pi + pa) and q(=|(pi — P2)) are, respectively, the total and relative 
momenta of the particle pair; and 0(q, r) is the relative two-particle wave function with 
r being their relative position, i.e., r = (r 2 — i^) — |(vi + v 2 )(t2 — ti)- This approach has 
been very useful in studying effects of nuclear equation of state and nucleon-nucleon cross 
sections on the reaction dynamics of intermediate energy heavy-ion collisions In the 
present paper, we use the Koonin-Pratt method to determine the nucleon-nucleon correlation 
functions in order to study the effect due to the momentum dependence of nuclear mean- 
field potential and the density dependence of nuclear symmetry energy on the spatial and 
temporal properties of nucleon emission source in intermediate energy heavy-ion collisions. 
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FIG. 6: (Color online) Emission rates of protons and neutrons as functions of time for different 
nucleon effective interactions. 

As an example, we study here central collisions of 52 Ca + 48 Ca at E = 80 MeV/nucleon. 
This particular reaction system with isospin asymmetry 5 = 0.2 can be studied at the future 
Rare Isotope Accelerator (RIA). In the present work, nucleons are considered as emitted 
when their local densities are less than po/8 and subsequent interactions do not cause their 
recapture into regions of higher density. In Fig. El we show the emission rates of protons and 
neutrons as functions of time for the SBKD, MDYI, and MDI interactions with soft and hard 
symmetry energies. It is clearly seen that there are two stages of nucleon emissions; an early 
fast emission and a subsequent slow emission. This is consistent with the long-lived nucleon 
emission source observed in previous BUU calculations j^J. For the momentum- independent 
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nuclear potential (SBKD), we find from Fig. El (a) that the hard symmetry energy enhances 
the emission of early high momentum protons (dash-dotted line) and neutrons (dotted line) 
but suppresses late slow emission compared with results from the soft symmetry energy 
(protons and neutrons are given by solid and dashed lines, respectively). The difference 
between the emission rates of protons and neutrons is, however, larger for the soft symmetry 
energy. Fig. El (b) shows results from the MDYI interaction which includes the momentum- 
dependent isoscalar potential but the momentum-independent symmetry potential. It is seen 
that the momentum dependence of isoscalar potential enhances significantly the nucleon 
emission rate due to the more repulsive momentum-dependent nuclear potential at high 
momenta. As a result, the relative effect due to the symmetry potential is reduced compared 
with the results shown in Fig. El (a). In Fig. El (c), we use the MDI interaction which 
includes momentum dependence in both isoscalar potential and symmetry potential. The 
momentum dependence of symmetry potential leads to a slightly faster nucleon emission but 
the symmetry potential effects are reduced. We note that the fraction of total number of 
emitted nucleons, i.e., before 200 fm/c in the IBUU simulations, is about 80% for the SBKD 
interaction but almost 100% for the MDYI and MDI interactions. The enhancement of 
nucleon emissions with momentum-dependent nuclear mean-field potential was also observed 
in previous calculations (37 1. It should be noted that the emitted nucleons in the present 
study are not exactly those observed experimentally since we have not included explicitly 



ragments in the IBUU simulations. 

34j], which takes into account final-state 



the production of clusters and intermediate mass 
Using the program Correlation After Burner 
nucleon-nucleon interactions, we have evaluated two-nucleon correlation functions from the 
emission function given by the IBUU model. Shown in Fig. [7| are two-nucleon correlation 
functions gated on the total momentum P of nucleon pairs from central collisions of 52 Ca + 
48 Ca at E = 80 MeV/nucleon by using the SBKD interaction with the soft and hard sym- 
metry potentials. The left and right panels are for P < 300 MeV/c and P > 500 MeV/c, 
respectively. Both neutron-neutron (upper panels) and neutron-proton (lower panels) corre- 
lation functions peak at q ~ MeV/c, while the proton-proton correlation function (middle 
panel) is peaked at about q = 20 MeV/c due to the strong final-state s-wave attraction. The 
latter is suppressed at q = as a result of Coulomb repulsion and anti-symmetrization of 
the two-proton wave function. These general features are consistent with those observed in 
experimental data from heavy-ion collisions j^. 
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FIG. 7: (Color online) Two-nucleon correlation functions gated on the total momentum P of 
nucleon pairs using the SBKD interaction with soft (filled squares) or stiff (open squares) symmetry 
energy. Left panels are for P < 300 MeV/c while right panels are for P > 500 MeV/c. 

For nucleon pairs with high total momentum, their correlation function is stronger for the 
hard symmetry energy than for the soft symmetry energy: about 24% and 9% for neutron- 
proton pairs and neutron-neutron pairs at low relative momentum q = 5 MeV/ c, respectively, 
and 21% for proton-proton pairs at q = 20 MeV/c. The neutron-proton correlation function 
thus exhibits the highest sensitivity to the density dependence in nuclear symmetry energy 
E sym (p). For nucleon pairs with low total momenta, the symmetry potential effects are weak. 
These results are again consistent with previous calculations |28l l29j]. 

What will happen to the symmetry energy effect on two-nucleon correlation functions 
if we include the momentum-dependent isoscalar potential in the IBUU model? This can 
bee seen from Fig. IHJ where they are shown for the same central collisions of 52 Ca+ 48 Ca 
at E = 80 MeV/nucleon by using the MDYI interaction with soft and hard symmetry 
potentials in the IBUU model. For nucleon pairs with low total momentum, shown in left 
panels, their correlation functions remain insensitive to the nuclear symmetry energy. For 
nucleon pairs with high total momentum, shown in right panels, their correlation function is 
again stronger for the hard symmetry energy than for the soft symmetry energy: about 17% 
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FIG. 8: (Color online) Same as in Fig. [7|but using the MDYI interaction with soft (filled squares) 
or stiff (open squares) symmetry energy. 

and 4% for neutron-proton pairs and neutron-neutron pairs at low relative momentum q = 5 
MeV/c, respectively, and 12% for proton-proton pairs at q = 20 MeV/c. Compared to the 
results from the SBKD interaction, the momentum dependence of isoscalar potential thus 
reduces the symmetry potential effects on two-nucleon correlation functions. This is due to 
the fact that the repulsive momentum- dependent potential enhances nucleon emissions and 
thus reduces the density effect on nucleon emissions, leading to a weaker symmetry potential 
effects on two-nucleon correlation functions. 

To see how the momentum dependence of nuclear symmetry energy affects the effect due 
to its density dependence on two nucleon correlation functions, we show in FigJUJthe results 
from the MDI interaction with the soft and hard symmetry potentials for central collisions 
of 52 Ca + 48 Ca at E = 80 MeV/nucleon. It is seen that the effect due to the stiffness 
of nuclear symmetry energy is again small for nucleon pairs with low total momentum as 
shown in left panels. For nucleon pairs with high total momentum, shown in right panels, 
their correlation function remains stronger for the hard symmetry energy than for the soft 
symmetry energy: about 11% for neutron-proton pairs at low relative momentum q = 5 
MeV/c and 9% for proton-proton pairs at q = 20 MeV/c. The symmetry potential effect on 
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FIG. 9: (Color online) Same as in Fig. 0but using the MDI interaction with soft (filled squares) 
or stiff (open squares) symmetry energy. 

the correlation function of neutron-neutron pairs is, however, very weak. Compared with 
results from the SBKD and MDYI interactions, the two-nucleon correlation functions from 
the MDI interaction are thus smaller, mainly due to the very small difference between its 
neutron and proton potentials, especially for higher momentum nucleons as shown in Fig. 
Hand Fig. H 



IV. THE T/ 3 HE RATIO 



Light cluster production has been extensively studied in experiments involving heavy-ion 
collisions at all energies, e.g., see Ref. j3| for a recent review. A popular model for descri bing 



the production of light clusters in these collisions is the coalescence model, ( ,l 
for a theoretical review, which has been used at both intermediate 



energies 



see Ref. 



6( 



and high 



7l\ . In this model, the probability for producing a cluster is determined by the 
overlap of its Wigner phase-space density with the nucleon phase-space distribution at freeze 
out. Explicitly, the multiplicity of a M-nucleon cluster in a heavy-ion collision is given by 
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N M = G / dY ix dq il ---dr iM _ 1 dq iM _ l { ^ p7( T h> • ' ' r *M-i> <J»j*-i))- ( 10 ) 

ii>i 2 >... >im 

In the above, r» i; • • • , r iM _ 1 and q^, • • • , qi Af _ 1 are, respectively, the M — 1 relative coor- 
dinates and momenta taken at equal time in the M-nucleon rest frame; pf is the Wigner 
phase-space density of the M-nucleon cluster; and (• • • ) denotes event averaging. The spin- 
isospin statistical factor for the cluster is given by G, and its value is 3/8 for deuteron and 
1/3 for triton or 3 He, with the latter including the possibility of coalescence of a deuteron 
with another nucleon to form a triton or 3 He (3]. Details about such calculation can be 
found in Ref. I3JJ. 
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FIG. 10: (Color online) Kinetic energy spectra in the center-of-mass system for deuteron (first 
row), triton (second row), and 3 He (third row) from central collisions of 52 Ca + Ca at E = 80 
MeV/nucleon by using the SBKD (first column), MDYI (second column) and MDI (third column) 
interactions with the soft (solid squares) or stiff (open squares) symmetry energy. 



Shown in Fig. El are the kinetic energy spectra in the center-of-mass system for deuterons 
(first row), tritons (second row), and 3 He's (third row) from central collisions of 52 Ca + 
48 Ca at E = 80 MeV/nucleon by using the SBKD (first column), MDYI (second column) 
and MDI (third column) interactions with the soft (solid squares) and stiff (open squares) 
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symmetry energies. For the SBKD interaction, the yields of these light clusters are very 
sensitive to the density dependence of nuclear symmetry energy, which is consistent with 
previous studies [^[y]]. Including momentum dependence in the isoscalar potential by using 
the MDYI interaction in the IBUU model reduces slightly the symmetry potential effect, 
although the yields of light clusters are increased significantly as a result of enhanced nucleon 
emissions from the momentum-dependent nuclear potential as shown in Fig. |H1 If we further 
include the momentum dependence in symmetry potential by using the MDI interaction, 
the final yields of light clusters become even less sensitive to the density dependence of 
symmetry potential than those from the MDYI interaction. This is due to the very small 
difference between the neutron and proton potentials and the momentum dependence of their 
symmetry potentials in the MDI interaction as discussed above. It should be noted that 
the coalescence model is expected to be a reasonable model for describing the production 
of light clusters with large kinetic energies as they are not contaminated by contributions 
from decays of heavy fragments, which are mainly of low kinetic energies and may not be 
negligible in intermediate energy heavy-ion collisions. 




(b) MDYI 

52 Ca + 48 Ca 

E=80 AMeV, b=0 fm 
■ Soft Sym. Pot. 
□ Hard Sym. Pot. 
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FIG. 11: (Color online) The t/ 3 He ratio as a function of cluster kinetic energy in the center-of-mass 
system for different interactions (a) SBKD, (b) MDYI, and (c) MDI with the soft (solid squares) 
and stiff (open squares) symmetry energies. The lines are drawn to guide the eyes. 



The isobaric yield ratio t/ 3 He is less model-dependent and also less affected by other 
effects, such as the feedback from heavy fragment evaporation and the feed-down from pro- 
duced excited triton and 3 He states. In Fig. ^2 we show the t/ 3 He ratio with statistical 
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errors as a function of cluster kinetic energy in the center-of-mass system for the SBKD, 
MDYI and MDI interactions with the soft (solid squares) and stiff (open squares) symmetry 
energies. For all nuclear potentials, it is seen that the ratio t/ 3 He obtained with different 
symmetry energies exhibits very different energy dependence. While the t/ 3 He ratio in- 
creases with kinetic energy for the soft symmetry energy, it decreases and/or increase weakly 
with kinetic energy for the stiff symmetry energy. For both soft and stiff symmetry energies, 
the ratio t/ 3 He is larger than the neutron to proton ratio of the whole reaction system, 
i.e., N/Z= 1.5. This is in agreement with results from both experiments and the statistical 
model simulations for other reaction systems and incident energies 

flBBBQ it is 

interesting to note that although the yield of light clusters is not so sensitive to the density 
dependence of symmetry potential for the MDI interaction, the t/ 3 He ratio shows very dif- 
ferent energy dependence for the soft and hard symmetry potentials. This is related to the 
different momentum dependence of symmetry potential in the MDI interaction, especially 
at low densities, as shown in Fig. HJ 



V. SUMMARY 



Through studying two-nucleon correlation functions and light cluster production using an 
isospin- and momentum- dependent transport model, we have investigated the effects due to 
the momentum dependence of isoscalar nuclear potential and also the symmetry potential 
on the space-time properties of nucleon emission source in heavy-ion collisions induced by 
neutron-rich nuclei at intermediate energies. It is found that the momentum dependence 
of both isoscalar nuclear potential and the symmetry potential influences significantly the 
space-time properties of nucleon emission source, leading thus to appreciable effects on two- 
nucleon correlation functions and light cluster production in these collisions. Specifically, 
the momentum dependence of nuclear potential reduces the sensitivity of two-nucleon cor- 
relation functions and the light cluster yield on the stiffness of nuclear symmetry energy. 
However, the t/ 3 He ratio is still found to be sensitive to the stiffness of symmetry energy 
after including the momentum dependence of nuclear potential. The study of these observ- 
ables in intermediate energy heavy-ion collisions thus offers the possibility to probe both 
the momentum dependence of nuclear symmetry potential and the density dependence of 
nuclear symmetry energy. 
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